International Journal of Electrical and Computer Engineering (IJECE) 
Vol. 12, No. 2, April 2022, pp. 1859~1868 
ISSN: 2088-8708, DOI: 10.1159 1/ijece.v12i2.pp1859-1868 O 1859 


A new approach for solving multi-pantograph type delay 


differential equations 


Nidal Anakira‘, Ali Jameel’, Mohmmad Hijazi*, Abedel-Karrem Alomari‘*, Noraziah Man? 
'Department of Mathematics, Faculty of Science and Information Technology, Irbid National University, Irbid, Jordan 
School of Quantitative Sciences, Universiti Utara Malaysia (UUM), Kedah, Sintok, Malaysia 
3Department of Mathematics, College of arts and Sciences Tabarjal, Jouf Univeristy, Sakakah, Saudi Arabia 
‘Department of Mathematics, Faculty of Science, Yarmouk University, Irbid, Jordan 


Article Info 


ABSTRACT 


Article history: 


Received Apr 14, 2021 
Revised Aug 16, 2021 
Accepted Sep 10, 2021 


Keywords: 


Laplace transform 

MRPSM 

Pade approximant 

RPSM 

System of multi-pantograph 
Type differential equation 


In this paper, a modified procedure based on the residual power series 
method (RPSM) was implemented to achieve approximate solution with 
high degree of accuracy for a system of multi-pantograph type delay 
differential equations (DDEs). This modified procedure is considered as a 
hybrid technique used to improve the curacy of the standard RPSM by 
combining the RPSM, Laplace transform and Pade approximant to be a 
powerful technique that can be solve the problems directly without large 
computational work, also even enlarge domain and leads to very accurate 
solutions or gives the exact solutions which is consider the best advantage of 
this technique. Some numerical applications are illustrated and numerical 
results are provided to prove the validity and the ability of this technique for 
this type of important differential equation that appears in different 
applications in engineering and control system. 


This is an open access article under the CC BY-SA license. 


Corresponding Author: 
Nidal Anakira 


Department of Mathematics, Faculty of Science and Technology, Irbid National University 


2600 Irbid, Jordan 


Email: alanaghreh_nedal@ yahoo.com 


1. INTRODUCTION 
In this study, the following system of multi-pantograph type DDEs will be considered. 


u',(t) = Bu, (t) + Alt Uy (440), u2 (@11t), ~, Un (@11t)) 
u'a (t) = P2u2 (t) + fo(t, uy (@21t), Up (@21t), «.., Un (@21t)) (1) 


Subject to the following initial conditions 
uj(to) = uj, to <t <T 


where, p;i, u; are finite constant, f;, are analytical functions such that 


0< Qij <1, Qij = 1,2, n. 


and u,(t), i = 1,2,..,n. are unknown functions that will be determined. In recent years, numerous 
numerical and approximate analytical studies related to the delay differential equations (DDEs) have been 
appeared especially the multi-pantograph type DDEs which is an important types of DDEs that appears in 
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various applications of engineering and branches of applied science such as electronic systems, population 
dynamics, dynamical systems, quantum mechanics. Obtaining analytical solutions of these kinds of DDEs are 
generally difficult. Therefore, in literature, many powerful and effective numerical procedures are used to 
approximate the solutions of multi-pantograph DDEs, for examples, optimal homotopy asymptotic method 
(OHAM) and its modifications are used to obtain approximate solutions for multi-pantograph time DDEs 
[1]-[7]. A domain decomposition method (ADM), differential transformation method (DTM) are also 
another approaches which they employed to solve multi-pantograph equations [8]. A collocation method is 
also one of the methods which is employed to find an approximate solution of a system of multi pantograph 
type delay differential [9]. Moreover, a number of solution methods have been proposed and employed for 
the approximate solutions of multi-pantograph type DDEs [10]-[13]. Mathematical modeling of physical- 
world phenomena leads to linear or nonlinear ordinary or partial differential equations in different branches 
of applied science and engineering. It is well-known that the solution of the nonlinear one is not easy either 
numerically or theoretically and even more difficult to establish a real model for the nonlinear problem. 
Therefore, a researcher has put their efforts and attention to develop and construct analytical and numerical 
methods to study and analyze the solution of the differential equations trying to provide highly accurate 
solution which is closed form to the analytic solution or give it. 

The residual power series method (RPSM) is one of an effective and easy power series procedure 
which was used in a large scale in the last few years for various forms of strongly linear and nonlinear 
differential equations and fractional differential equations without any limitations or restrictions such as 
linearization, discretization and perturbation [14]—[20]. Also, it has been successfully employed to handle the 
numerical solution of a highly nonlinear singular Lane-Emden equation type [21]. The fundamental 
motivation of this paper is based on improving the accuracy of the RPSM using an alternative procedure for 
the first time that modify the series solution of the RPSM throughout using the Laplace transformation to the 
truncated RPSM solution then convert the transformed series into a meromorphic function by using Pade 
approximants, lastly, we apply the inverse of the Laplace transformation to get the required solution of the 
given problem. This method is easy, and does not requires big efforts to achieve accurate results with high 
performance. 

The rest of this paper is arranged as follows. Section 2 is conducted to illustrate the fundamental 
concept of the RPSM along with the method analysis besides to preliminary of the Pade approximants. In 
section 3 numerical examples are given to prove the capacity of the discussed procedure by comparing the 
results which indicate that only a few terms are required to deduce the exact solutions. Finally, conclusions of 
this work are formulated in the last section. 


2. RESEARCH METHOD 
2.1. Residual power series method (RPSM) 

The aim of this part is to give an overview and some properties of the RPSM procedure, which will 
be employed in order to get series solution of a system of multi-pantograph type DDEs [22]—[24]. The RPSM 
provide us the solutions of the given system in a form of power series expansion about the initial point t = to 
to receive our objective, we assume the solutions in the form of (3): 


uilt) = Dn=o Uim(t), i = 1,2, n, (3) 


where uUj.(t) are the order of approximations such that Uim(t) = Cim(t — to)”. Then, by substituting the 
initial guesses Uj 9(to) = u™ (to) = Cio Which are the initial conditions (2), using m = 0, into 
u;(t),i = 1,2, yield to the approximate solution for the system of multi-pantograph equations 
uj(t) = ui olto) + Ym=o Uim(t), i = 1,2, ...n, whereas ujm(t) for m = 1,2,...,k, can be obtained from the 
following k"-truncated series. 


Uj R(t) = Ds Cim (t = to)”, i = 1,2, “N, (4) 
To use the RPSM, we write the system of (1) and (2) in the form given: 


u'i (t) — piui Çt) — fi uy (airt), U2 (izt), «Un (@int)) 
u'a(t) — Bout) — fo(t, ur (@it), U2 @izt), ..., Un(@int)) (5) 


The k" and the oo" residual functions are given respectively, in the following manner 
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Resi (t) = u'ik (t) — Buin (t) — fit, u1 k (ait), Uzlat), Unk (int), i = 1,2,...,n (6) 
and 
Res? (t) = lim Resi (t) =u':(t) — piut) - 
fit, Uy (@j1t), Uz (ait), Un (jn t), i = 1,2,..,n (7) 


Clearly, we see that Res;°(t) = 0 for each t E (to,T), are infinitely differentiable functions at t = tọ. 
Moreover, a Res? (to) = T Te =0,m= i 2,...,k, which is a fundamental rule in the RPSM 


and its applications. In particular, © a L Res? (to) = ~ Resk (to) = 0,i = 1,2,...,n,k =1,2,.... 


Then by substituting the k"truncated series Ui p(t) into (6) we get 


Io 


Res (t) = Emsi MCi m(t — to) ™1 
—B; no Cim(t — to)” — filt, Lmao cım (ait — to)”, 
Eh-o cM (ant — tp) iss Dad Cnm (Aint — ty) i = Dd (8) 


Hence, depend on Res} (ty) = 0, i = 1,2, ...,n. Using t = tọ = 0 and k = 1 into (8) yields to 
Cia = PiCio + fi(to, C10, C20 ~ Cn,0) = Bio + fi(0, Tio) i =1,2,...,n (9) 


where fi(0, Ui) = fi(0, u0, U20» s Uno): Using (4), the first approximate solutions of the system (1) and 
(2) will be written in the following form (10). 


uia (Ë) = ttio + (Bitio + fi(0, Tio) ) (t — to), i = 12,40 (10) 


The second-order approximate solutions can be obtained by using k=2 and tọ=0 such that 
Uj2(t) = X2,-0 Cimt™. Then, by differentiate with respect to t, we have (11). 


2 


d d 
= Res? (0) = 2ci2 — Bica - E > amalit", 
m=0 
Dao maat beatae DES Lean (11) 


Based on the fact that the S ~ Resk (0) = 0,i = 1,2, ...,n, the values of c; z are given by (12), 


2 


1 d 
Gia = 5 Picat Lilt, > mare”, 
m=0 


Ce s Caine leo = 1,2, n (12) 


Therefore, the second approximate solutions will be (13), 


Uiz = Uio + (Biuio + fi(0, Tio) ) (t — to) + 
5 (Bicia + 91(0, T4 )Ct — to), i = 1,2, n (13) 


where g9i(0, Uj a= Tine, ae o C1, m4 M Ema o C2 maint™, Dorp knee Cnmaint” le=0) i = 1, 2, ay NM and 
Ci, are given in (9) fot i=1,2,...,n. We Sola the same procedure to generate a sequence of approximate 
solutions u;,(t) for the system (1) and (2) Furthermore, high accuracy can be achieved by increasing the 
order of the approximations. For this regards, we let Rem*(t) to denote the difference between u;(t)and its 


kth Taylor polynomial, that is, here the functions Rem¥ (t) represent the kth remainder of the Taylor series 
of ui (t) 


Remi (t) = uj(t) + uilt) = Leiner UT (to) (t — to)” „i = 1,2, nn (14) 


et 
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2.2. Padè approximation 
The [L/M] Padé approximants [25]-[27] of a function u(x) is defined by: 


P 3 mG 


where P,(t) and Qy(t) are polynomials of degrees at most L and M, respectively. The general form of the 
power series is 


u(t) = Xiz at’. 


The coefficients of the P, (t) and Qu (t) polynomials are obtained from: 


_ Pr) _ L+M+1 
u(t) oT O(t ) (15) 


O is multiplying by a nonzero constant the fractional 
M 


values remain unchanged, then we can define the normalization condition as (16): 


when the fraction of the numerator and denominator 


Qu(0)=1 (16) 


It can be noted that P, (t) and Qy(t) have no public factors. If we repsents the coefficient of P, (t) and Qy(t) 
as (17): 


l P(t) = po + pit + pat? + © + pit” 


(17) 
Qu(t) = qo + qıt + q2t? + © + qut™ 


Then, by (16) and (17), we can multiply (15) by Qu (t), to linearizes the coefficient equations. We can write 
out (18) in more detail as (18): 


at+ı + aLqı t+ + aL-m+1qm = 0 
al+2 + aL+1q1 +7 + aL-m+2qm = 0 


apim + aL+Mm-1q1 +’ + aqm = O (18) 


To obtain the solutions of these equations, we begin with (18), which is consider as a set of linear 
equations for all of the unknown q's. Once the q's are known, then (19) gives an explicit formula for the 
unknown p's, which complete the solution. If (18) and (19) are non-singular, then we may solve them direct 
and we obtain (21) [23], where (21) holds, and if the lower index on a sum exceeds the upper, the sum is 
replaced by zero: 


ao = Po 


ao + adi = Py 
az + a1qı + aAgq2 = P2 


ay + ap-141 +*+ aoqL = PL (19) 


aL-M+1 aL-M+2 we aL+1 
det . 
aL aL+1 os aL+M 
L j L j L j 
L] Èj=m aj—mX’ Èj=m-1 âj-m+1X' dito ajx? 
M aL-M+1 aL-M+2 nee aL+1 
det ; 
ay alti ua aL+M 
xM xM-1 1 (20) 
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Therefore, the Pade approximants diagonal matrix of different order can be obtained using Mathematica 


software or like MATLAB and son. 


3. RESULTS AND DISCUSSION 
3.1. Example 1 
The following system of multi-pantograph type delay differential equations considered [21]. 


w — u(t) + u(t) — u (5) = g0, 


u(t) + u (O) + u2(t) + uz ($) = 8:0 (21) 

u,(0)=1 u,(0)=1 (22) 
where 

gE ete  g)(t)=e'+e2 23) 


To solve this system by using RPSM, we start with the initial conditions u4 o = 1 and uzo = 1, as 
an initial guesses, then, we consider the kth-truncated series solutions of u4 g(t) and uz g(t) in form (24): 


welt) = Eh-0 Gmt” = 1+ cat +t? + cgt? He,t", 
uz p(t) = ee Comt™ = 1 + Coit + C2 2t? + C23 t? +e Cart" (24) 


Now, the unknown coefficients cj m,i = 1,2. m = 1,2,...,k can be obtained by formulated the following k* 
residual functions 


Rest (t) = Frizi MC mt”! mz er Cimt” + ino Comt™ 
t _ r 

5 Dnd mG” —e™ ee, 

Res} (t) = Xh=1 MC mt”! + Snö C1mt™ + Daci Com t™ 


t =E 
rE =e! hes (25) 


The first-order approximation solution u,,(¢) = 1 — t and u, (t) = t can be obtained based on the 
above residual function by using k=1 to get 


3 t 
Resi (t) = c11 + (cos — Zer) t-e'+e2-1 


É 
Resį (t) = c2, + (cia + ea) t-e2-e'+3 (26) 
Based on the fact that £ Res} (0) = £ Res}(0) = 0, we have c,, = 1 and c,, = —1, which yields to the 
first order approximate solution u,, = 1 + t and uz, = 1 — t. The values of the coefficients c, and c2 2 can 
be found by differentiate both sides of (25) with respect to t when k=2 and then substituting t=0 to be c12 = 


1 1 : ; : : 1 
5 and c22 = — 3 10 obtain the following second order approximate solution u12 = 1+¢t+ z6 and uz, = 1+ 


1 ; ; . 
t=5t By follow the same proceeder, the 5" order approximate solutions for the given system become 


$3 tt. #255 8 
-+-+ 
3l 4! 5! 6 

th tt E ENE £8 


3! 4! 5! 6! 


2 
Ue) =1+t+2+ 


2 
u(t) =1-t+ a (27) 
which leads to the analytic form of the solution in the limit of infinity terms of the order of the 
approximation. To improve the efficiency of the solutions procedure and get high accuracy we will use an 
effective and powerful modification based on RPSM truncated series solutions by employing the Laplace 
transform to the first four terms of (27) to obtain. 
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Pig benr krl ces SOE SE: T 
Lus O}= stata ta LsO ataa (28) 


s s2 s 
Consider s = E, yields 
Lfu,s(t)} =z7+2z? +z 427%, Lfu s (t)} =z—z*+z3 —z* (29) 


Then, using the Pad approximants of Él. 


Lu = — (30) 
Noting that z = *, and by using the inverse of the Laplace transform, we have 


u(t) = ef 
u(t) = e™t (31) 


which is the exact form of the above system. 


3.2. Example 2 
Consider the following nonlinear system of multi-pantograph type delay equation [21]. 


ui (t) = u (t) — e™% cos(Ż)uz © ei £) sin £), 
u) = e? (£) - u3 È) ae 


2 


Subject to given initial conditions 

u,(0)=1, u,(0) =0 (33) 
To obtain the approximate solution using RPSM, we will start with the initial conditions 

Uzo(t) = 1, Ugo(t) =0 (34) 
The kth-truncated series formula (3) for this example by using (17) is 


ui k(t) — aa Cimt™ = 1 + Cuil + Cy 2t? + czt’ +e , Cirt”, 
uz x(t) = er Comt” e 1 + Coit + C2 2t? + Co3t? +e. „Capt“ (35) 


whereas the k" residual functionRes* (t), where i = 1,2, is 


a _ t t 
Resp) = Xk,- Mei, t+ eb Cmt” +e™t cos(>) hat Com G)” 
3 


+2e(-3) cos È) sin È) Desg Cum 9”, 
ReO = Sh aimee E O e TS Oe) (36) 


when k = 1, we get the first order RPS approximate solutionu,,(t) = 1 — t and u,,(t) = t. The 2nd 
orderRPS approximate solution can be obtained by differentiate (36) and using t = 0,c1o =1 and 


a (Res}*(0) = 0); 


C29 = 0. We follow the same procedure up to the k" order based on the ee 


i = 1,2,10th order approximate solution is 


tt 2.6? t7 t8 
u t)=1-t — = 
110(t) F 3 6 t 30 630 + 2520" 
t3 t5 t7 t? 37 
u t)=t+—- + = 
110(t) 6 120 5040 362880 (37) 
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which leads to the analytic form of the solution in the limit of infinity terms of the order of the 
approximation. To improve the efficiency of the solutions procedure and get high accuracy we will use an 
effective and powerful modification based on RPSM truncated series solutions by employing the Laplace 
transform to the first four terms of (37) to obtain. 


1 1 r 4 4 
L{u,s(@} E eo ee F st 53 + 56 
1 1 1 1 1 
LwsO}=s-ata-ata (38) 


For the simplicity purpose, let s = 5, yields to 


L{u,s()} = z — z? + 224 — 42° + 42° 
L{uy5(t)} = z? + z4 — 26 + z8 + 219 (39) 


Now, we use the Pad approximants of ÉI. 


2 2 
Z°+Z w= Z (40) 


~ 22242241’ z241 


Lu 


; 1 : ; ; 
Consider z = > and by using the inverse of the Laplace transform to obtain 


u (t) = e™ cos(t) 
u,(t) = sint, (41) 


which is the exact form of the above system. 


3.3. Example 3 
Consider the nonlinear system of multi-pantograph type delay equation [21]. 


u,'(t) = 2u, £) + uz (t)tcos (<), 
u,'(t) = 1 — tsint — 2u2 (<) 
uz (t) = u2(t) — u,(t) — tcost (42) 


Subject to the given initial conditions. 
u,(0) = —1 u,(0) = 0 u3(0) = 0 (43) 


To obtain the approximate solution using RPSM, we start with the initial approximation 
Uz,9(0) = —1 uz 9(0) =0 uzo(0) = 0. Then we follow the same procedure which applied in the previous 
examples, to obtain the following truncated RPS series approximation 


t* te t8 t10 


2 

— eR m — t 

u,,(t) = =0 limt” = -14+-—--—+ — + 
1K(t) = Ym=0 Cim 2 24 720 40320 3628800 


3 5 7 9 
Z Pia t t t t 
Uz KCC) = Ym=0 C2mt™ = t — 2 24 720 40320 ” 
t3 t5 t7 9 


Uz k(t) = Lin=0C3mt™ =t +s, (44) 


6 | 120 5040 | 362880 
which leads to the analytic form of the solution in the limit of infinity terms of the order of the 
approximation. To improve the efficiency of the solutions procedure and get high accuracy we will use an 
effective and powerful modification based on RPSM truncated series solutions by employing the Laplace 
transform to the first four terms of (44) to obtain: 


ae Oe ee ee 
L{uy4(t)} = Saty oh big at 


9 7 5 3 1 
L{uz 0} =a ath at 
1 1 1 1 1. 
L{uz4(t)} = mio at 2 atz (45) 


A new approach for solving multi-pantograph type delay differential equations (Nidal Anakira) 


1866 O ISSN: 2088-8708 
For the purpose of simplicity, let= *, yields. 

L{u,5(t)} = -2° +27 - z" +23 —z, 

L{uz5(t)} = 921° — 7z8 + 52° — 324 + 2?, 


L{uz s (t)} = 21° — z8 + z6 — 24 + 2? (46) 


Now, we use the Pad approximants of BE 


Z 


L, =— 
u Z241 

n z2—z4 
U` 74-22241’ 

ly = = 47 
uu 2241 ( ) 


In (47) recalling that z = Z, and then applying the inverse Laplace transform to obtain 


u,(t) = —cos(t) 
u(t) = tcos(t) 
u3(t) = sin(t) (48) 


which is the exact form of the above system. From the numerical results obtained, we saw that the proposed 
procedure provides us high accurate solutions closed to the exact one. This can be obtained by using a few 
numbers of iterations of the standard RPSM solutions, this advantage overcomes the difficulties that its need 
to compute more iterations to get more accurate solutions to improve efficiency of the standard RPSM that is 
noted from the absolute errors represented by Figure 1 to Figure 3. 


0.0001 0.00008 
0.00006 
ral 
0.00004 


0.00002 


(a) (b) 


7. 10° 
6. 10° 
0.000015 | 
5.10 
4.10° 
>, 0.00001 os 
3.10 
2. 10° 
5. 10° 

1. 10° 

es ot: 

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

x x 
(a) (b) 


Figure 2. Absolute errors related to (a) u4 (t) and (b) uz (t) of example 2 
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4. 


5 


4 


3 


1” rma 8.108 


. 10 8 ; 
6. 10° 


. 10" 


1019 


2.108 
. 10 19 


oo 02 04 ~ 06 08 L0 wo x 04 06 os o 


0.0 0.2 0.4 0.6 0.8 1.0 


Figure 3. Absolute errors related to (a) u4 (t) and (b) u3 (t) of example 1, and (c) u3 (t) of example 3 


CONCLUSION 
In this research study, a new procedure based on the RPSM was proposed to solve system of 


multi-pantograph type DDEs. This procedure is effective, reliable and has a distinct advantage over the other 
methods represented by it is valid for the strongly nonlinear problems. It has been shown throughout the 
illustration examples and the comparison with of the numerical results reported in the literature that the 
MRPSM procedure has the sufficient to obtain the exact analytical solution using only few terms of the 
RPSM truncated series solution. This leads to conclude that this procedure is a powerful approach and a 
promising tool for solving this kind of differential equations. 
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